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1 Introduction 

We present a study and differential cross section measurements of photons emitted in decays 
of Z bosons produced at a hadron collider. Such radiative decays of the Z boson were noted 
in the very first Z boson publications of UA1 and UA2 dEIl, but subsequently have not been 
given a detailed study in hadron colliders. In 2011, the CERN LHC delivered pp collisions 
at i/s = 7TeV, and data corresponding to an integrated luminosity of 4.7 fb 1 were collected 
with the CMS detector. From these data, we select a sample of events in which a Z boson 
decays to a pair and an energetic photon. We measure the differential cross sections 

da/ dEj with respect to the photon transverse energy Ej and da/dAR /l7 with respect to the 

separation of the photon from the nearest muon. Here, AR /l7 = sj (c/y — cp 7 ) 2 + — ?/ 7 ) 2 , 

where (p is the azimuthal angle (in radians) around the beam axis and ?/ is the pseudorapidity. 
The cross sections include contributions from the Z resonance, virtual photon exchange, and 
their interference, collectively referred to as Drell-Yan (DY) production. 

Photons emitted in Z boson decays, which we call final state radiation (FSR) photons, can be 
energetic (tens of GeV) and well separated from the leptons (by more than a radian). Quantum 
electrodynamics (QED) corrections that describe FSR production are well understood. Quan¬ 
tum chromodynamics (QCD) corrections modify the kinematic distributions of the Z boson; in 
particular, the Z boson aquires a nonzero component of momentum transverse to the beam: 
cjj > 0. The FEWZ program calculates both QCD and QED corrections for the DY process 0. 
However, it does not include mixed QCD and QED corrections; the required two-loop inte¬ 
grals are technically very challenging, and progress has been made only recently JU. In prac¬ 
tice, event generators employing matrix element calculations matched to parton showers must 
be used 10 16). It is the goal of this analysis to establish the quality of the modeling of FSR 
over a wide kinematic and angular range. The results will support future measurements of 
the W mass, the study of Z +7 production, and searches for new particles in final states with 
photons. 

In an attempt to compare photons emitted close to a muon (a process that is modeled primar¬ 
ily by a partonic photon shower) and far from the muons (which requires a matrix element 
calculation), we measure da/ dEj for 0.05 < AR „ 7 < 0.5 and 0.5 < A R in < 3. Furthermore, 
since the size of the QCD corrections varies with the transverse momentum of the Z boson, 
we measure da/ dEj and da/ dA R in for qj < 10 GeV and qj > 50 GeV, where qj is defined 
as the vector sum of the transverse momenta of the two muons and the photon. These cross 
sections are defined with respect to the fiducial and kinematic requirements detailed below; 
no acceptance corrections are applied. Nonetheless, we do correct for detector resolution and 
efficiencies. 

This article is structured as follows. We briefly describe the CMS detector and the event samples 
we use in Section [2j The details of the event selection are given in Section [3] Background 
estimation and the way we unfold the data distributions are discussed in Sections [4] and [5] 
We discuss the systematic uncertainties in Section [ 6 ] and report our results and summarize the 
work in Sections [7] and [ 8 ] 


2 The CMS detector and event samples 

A full description of the CMS detector can be found in Ref. 13; here we briefly describe the 
components most important for this analysis. The central feature of the CMS experiment is 
a superconducting solenoid that provides an axial magnetic field of 3.8 T. A tracking system 
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3 Event selection 


composed of a silicon pixel detector and a silicon strip detector is installed around the beam 
line, and provides measurements of the trajectory of charged particles for \t]\ < 2.5. After 
passing through the tracker, particles strike the crystal electromagnetic calorimeter (ECAL) 
followed by the brass and scintillator hadron calorimeter. The solenoid coil encloses the tracker 
and the calorimetry. Four stations of muon detectors measure the trajectories of muons that 
pass through the tracker and the calorimeters for \p\ < 2.4. Three detector technologies are 
employed in the muon system: drift tubes for central rapidities, cathode strip detectors for 
the forward rapidities, and resistive-plate chambers for all rapidities. Combining information 
from the muon detectors and the tracker, the transverse momentum (pj) resolution for muons 
used in this analysis varies from 1 to 6%, depending on ?/ and pj |8|. The Ej of photons and 
electrons is measured using energy deposited in the ECAL, which consists of nearly 76 000 lead 
tungstate crystals distributed in the barrel region (|;/1 < 1.479) and two endcap regions (1.479 < 

| ?/1 < 3.0). The photon energy resolution is better than 5% for the range of Ej pertinent to this 
analysis |9j- Events are selected using a two-level trigger system. The level-1 trigger, composed 
of custom-designed processing hardware, selects events of interest based on information from 
the muon detectors and calorimeters HlQll . The high-level trigger is software-based, running 
a simpler and therefore faster version of the offline reconstruction code on the full detector 
information, including the tracker lUTI . 

Simulated data samples are used to design and verify the principles of the analysis. They are 
also used to assess efficiencies, resolution, and backgrounds. The signal process is simulated 
using the POWHEG (VI.0) lH2l event generator with PYTHIA (V6.4.24) llT3l used to simulate 
parton showers and final state photons (referred to in what follows as POWHEG+PYTHIA). This 
combination is also used for tt and diboson (WW, WZ, ZZ) production. The CT10 Ifl4l parton 
distribution functions are used. The Z2 parameter set jT5l is used to model the underlying 
event in PYTHIA, and the effects of additional pp collisions that produce signals registered 
together with the main interaction are included in the simulation. 

The response of the detector is simulated using Geant4 [16J. The simulated events are pro¬ 
cessed using the same version of the offline reconstruction code used for the data. 


3 Event selection 

The data are recorded using a trigger that requires two muons. One muon is required to have 
pj > 13 GeV, and the other, pj > 8 GeV. This trigger has no requirement on the isolation of the 
muons. 

Events with a pair of oppositely charged, well-reconstructed, and isolated muons and an iso¬ 
lated photon are selected. The kinematic and fiducial requirements for selecting events are 
based wholly on the muon and photon kinematic quantities, and are summarized in Table [T| 
As explained below, we use the dimuon mass M w to define a "signal region" that is rich in 
FSR photons, and a "control region" that is dominated by background sources of photons. 

Muons are selected in the manner developed for the measurements of the DY cross section lH7| . 
They must be reconstructed using an algorithm that finds a track segment in the muon detec¬ 
tors and links it with a track in the silicon tracker, and also through an algorithm that extrap¬ 
olates a track in the silicon tracker outward and matches it with hits registered in the muon 
detectors. We select the two highest pj muons (which we will call "leading" and "trailing"), 
and ignore any additional muons. These two muons must have opposite charge. The leading 
muon must satisfy the requirements pj > 31 GeV and // < 2.4, while the trailing muon must 
satisfy pj > 9 GeV and \t]\ < 2.4 to ensure good reconstruction efficiency. A vertex fit is per- 
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formed to the two muon tracks, and the x 2 probability of the fit must be at least 0.02. We define 
the difference between n and the opening angle of the two muons as the acollinearity a, and 
remove a very small region of phase space where a. is less than 5 mrad to reduce contamination 
by cosmic rays to a negligible level. 

Photons are reconstructed using the particle flow (PF) algorithm lfT8lfT9l that uses clustered en¬ 
ergy deposits in ECAL. The PF algorithm allows us to reconstruct photons at relatively low Ej 
and to maintain coherence with the calculation of the isolation observables described below. 
Photons that convert to electron-positron pairs are included in this reconstruction. Events se¬ 
lected for this analysis must have at least one photon with Ej > 5 GeV, and the separation of 
this photon with respect to the closest muon must satisfy 0.05 < AR /(7 < 3. Studies using the 
simulation show that photons with AR W7 < 0.05 are difficult to reconstruct reliably due to the 
energy deposition left by the muon, and no signal photons with A R }n > 3 are expected. If 
an event has more than one photon satisfying this selection criteria, we select the one with the 
highest Ej. In events in which one photon is selected, a second photon is present 15% of the 
time; however, these extra photons are expected to be mostly background, since the fraction 
of events with a second FSR photon in simulation is approximately 0.5%. More details about 
these background photons are given in Section |4] 

Table 1: Summary of kinematic and fiducial event requirements 


Object 

Requirement 

heading muon 

pj > 31 GeV and \rj\ < 2.4 

Trailing muon 

pj > 9 GeV and \tj\ < 2.4 

Acollinearity 

a. > 0.005 radians 

Photon 

£ t > 5 GeV, 

\f]\ < 2.4 but not 1.4 < \t]\ < 1.6; 
0.05 < A R in < 3 

Signal region 

30 < M m < 87 GeV 

Control region 

89 < M m < 100 GeV 


All three particles emitted in the Z boson decay—the two muons and the photon—are usually 
isolated from other particles produced in the same bunch crossing. We can reduce backgrounds 
substantially by imposing appropriate isolation requirements. The isolation quantities, /„ for 
the muons and 1 7 for the photon, are based on the scalar pj sums of reconstructed PF particles 
within a cone around the muon or photon direction. The cone size for both muons and photons 
is A R = 0.3. The muon pj is not included in the sum for / ;i , and the photon Ej is not included in 
the sum for 1 7 ; these isolation quantities are meant to represent the energy carried by particles 
originating from the main primary vertex close to the given muon or photon. For a well- 
isolated muon or photon, / ;i or I 7 should be small. 

Special care is taken to avoid inefficiencies and biases occurring when the FSR photon falls 
close to the muon; in such cases the muon and the photon may appear, superficially, to be 
nonisolated. To avoid this effect, we exclude any PF photon from the muon isolation sum. Fur¬ 
thermore, since the photon can convert and produce charged particle tracks that cannot always 
be unambiguously identified as an e ! e pair, we exclude from the muon isolation sum charged 
tracks that lack hits in the pixel detector or that have pj < 0.5 GeV. Finally, any particle that 
lands in a cone of A R < 0.2 around a PF photon is excluded from the muon isolation sum. After 
these modifications to the muon isolation variable, the efficiency of the isolation requirement is 
flat (98%) for all AR ; , 7 and is higher than the efficiency of the unmodified isolation requirement 
by about 0.5%. Adding these modifications does not significantly increase the backgrounds. 
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4 Background estimation 


The instantaneous luminosity of the LHC was sufficiently high that each bunch crossing re¬ 
sulted in multiple pp collisions (8.2 on average). The extraneous pp collisions are referred to 
as "pileup" and must be taken into consideration when defining and calculating the muon and 
photon isolation variables. Charged hadrons, electrons, and muons coming from pileup can 
be identified by checking their point of origin along the beam line, which will typically not 
coincide with the primary vertex from which the muons originate. When summing the contri¬ 
butions of charged hadrons, electrons, and muons to the isolation variable, those coming from 
pileup are excluded. This distinction is not possible for photons and neutral hadrons, however. 
Instead, an estimate I p of the contribution of photons and neutral hadrons to the sum is made: 
we use one-half of the (already excluded) contribution from charged hadrons within the isola¬ 
tion cone. This estimate is subtracted from the sum of contributions from photons and neutral 
hadrons; if the result is negative, we then use a net contribution of zero. 

We designate by 4± the sum of pj for charged particles that are not excluded from the isolation 
sum. We let I em and / h o stand for the sums over the pj of all photons and neutral hadrons, and 
Ip for the estimate of the pileup contribution to I em and I h o. The muon isolation variable is, 
then, 

4 = (4± + 4°) / Pt- (1) 

Note that the sum is normalized to the pj of the muon. We require 4 < 0.2 for both muons. 

The photon isolation variable is calculated as above, except that the muons are not included in 
the sum, and there is no special exclusion of charged tracks near the photon: 

f 7 = 4± + max(I em + 4o - Ip, 0). (2) 

We require I 7 < 6 GeV. 

The emission of FSR photons in Z boson decays reduces the momenta of the muons. Con¬ 
sequently, the dimuon mass M (/// tends to be lower than Mz, the nominal mass of the Z bo¬ 
son. Simulations indicate that, for most of the signal, M,,,, < 87 GeV, due to the requirement 
Ej > 5 GeV for the photon. They also show that the M //; , distribution for radiative decays 
Z —> ends at ~ 30 GeV. A requirement > 30 GeV also helps to avoid a 

kinematic region in which the acceptance is difficult to model. Therefore, our signal region is 
defined by 30 < M ((/i < 87 GeV. We also define a control region by 89 < M ;i/ , < 100 GeV, 
where the contribution of FSR photons is quite small (below 0.5%). The numbers of events we 
select are 56 005 in the signal region and 45 277 in the control region. 

4 Background estimation 

Nearly all selected events have two prompt muons from the DY process. Backgrounds come 
mainly from "nonprompt" photons, which may be genuine photons produced in the decays of 
light mesons (such as tt° and ;/), a pair of overlapping photons that cannot be distinguished 
from a single photon, and photons from pileup. We study these backgrounds with simulated 
DY events and apply corrections so that the simulation reproduces the data distributions, as 
described in detail below. 

Some events come from processes other than DY, such as tt and diboson production. These 
backgrounds are very small and are estimated using the simulation. Similarly, a small back¬ 
ground from the DY production of X 1 x is also estimated from simulation. The background 
from multijet events, including events with a W^ —» }i ± v decay, is estimated using events with 
same-sign muons. Backgrounds from simultaneous nonprompt muon and nonprompt photon 
sources are negligible. The composition of the signal sample is given in Table [2] 
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Table 2: Composition of the signal sample. The simulation has been tuned to reproduce the 
data in the control region. 


Process 

Fraction 

Signal 

77.1% 

DY with a nonprompt photon 

9.5% 

Pileup 

11 .2% 

tt 

0 .6% 

t + t - 

0.3% 

Dibosons 

0 .2% 

Multijets 

1 .1% 


The control region is dominated by nonprompt photons whose kinematic distributions (Ej, 
if, A/% 7 ) are nearly identical to nonprompt photons in the signal region. Quantitative com¬ 
parisons of data and simulation revealed significant discrepancies in the control region that 
prompted corrections to the simulation, which we now explain. 

The powheg+pythia sample does not reproduce the number of jets per event well, so we 
apply weights to the simulated events as a function of the number of reconstructed jets in each 
event. For this purpose, jets are reconstructed from PF objects using the anti-fcj algorithm |20| 
with a size parameter R = 0.5. We consider jets with pj > 20GeV and \rj\ < 2.4 that do not 
overlap with the muons or the photon. 

Studies of f 7 for events in the control region reveal small discrepancies in the multiplicity and 
pj spectra of charged hadrons included in the sum. We apply weights to the simulated events 
to bring the multiplicities into agreement, and we impose pj > 0.5 GeV on charged hadrons. 
The simulated I 7 distributions match those in data very well after applying these weights. 

Finally, the Ej and if distributions of nonprompt photons in the simulation deviate from those 
in the data. We fit simple analytic functions to the ratios of data to simulated Ej and if distri¬ 
butions and define a weight as the product of those functions. We check that this factorization 
is valid (i.e., that the Ej correction is the same for different narrow ranges of if, and vice versa). 

After these three corrections (for jet multiplicity, for the spectrum of charged hadrons in the 
isolation sums, and for the Ej and if of the nonprompt photon), the simulation matches the data 
in all kinematic distributions in the control region, an example of which is shown in Fig. [Tj. left. 
The total change in the background estimate due to these corrections is approximately 5 to 10%. 
We use the simulation with these weights to model the small background in the signal region 

(Fig-0] ri S ht )- 

Given the definition of the signal region, the contribution of photons emitted in the initial state 
is very small (on the order of 4 x 10~ 4 as determined from the POWHEG+PYTHIA sample) and 
is counted as signal. 


5 Correcting for detector effects 

Our goal is to measure differential cross sections in a form that is optimal for testing FSR cal¬ 
culations. Therefore we are obliged to remove the effects of detector resolution and efficiency 
(including reconstruction, isolation, and trigger efficiency). The corrections for the muons fol¬ 
low the techniques developed for the DY cross section measurements fll7l . The corrections for 
photons are applied using an unfolding technique, as discussed in this section. 

We apply small corrections to the muon momentum scale as a function of muon pj, if ?( , and 
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6 Systematic uncertainties 


4.7 fb 1 (7 TeV) 



Photon Isolation (I) [GeV] 


4.7 fb -1 (7 TeV 



Photon Isolation (I) [GeV] 


Figure 1: Distributions of the photon isolation variable f 7 for the control region (left) and for 
the signal region (right) after all corrections have been applied. The bottom panels display the 
ratio of data to the MC expectation. The requirement for FSR photons is I 7 < 6 GeV. 

(pf, J2TJ; they have almost no impact on our measurements. The muon reconstruction efficiency, 
(taken from simulation and corrected to match the data as a function of pj and ?/,,) is taken into 
account by applying weights on a per-event basis. We do not correct for the approximately 
0.5% increase in the isolation efficiency coming from the way we handle FSR photons falling 
within the muon isolation cone. 

The energy scale and efficiencies for photons are more central to our task. Most PF photon 
energies correspond to the true energies within a few percent. However, in about 13% of the 
cases the photon energy is significantly underestimated. The simulation reproduces this effect 
very well. We construct a "response" matrix that relates the PF energy to the true energy as a 
function of ij 7 and AR^. The angular quantities ;/ 7 and AR f , 7 are themselves well measured. 
We use the iterative D'Agostini method of unfolding Il22l as implemented in the ROOUNFOLD 
package |[23l . By default, we unfold in the three quantities Ej, )/ 7 , and A R flJ simultaneously 
after subtracting backgrounds; as a cross check we also unfold the Ej and AK„ 7 distributions 
one at a time, and we also use a single-value decomposition method |[24ll . All results are con¬ 
sistent with each other. To verify the independence of the unfolded result on the assumed 
spectra, we distort the FSR model in an arbitrary manner when reconstructing the response 
matrix. The unfolded result is no different than the original one we obtained. A closure test in 
which the simulation is treated as data and undergoes the same unfolding procedure indicates 
no deviation greater than 1.5%. 

The unfolding procedure corrects for the photon reconstruction and isolation efficiency. It also 
corrects for bias in the PF photon energy assuming that such a bias is reproduced in the sim¬ 
ulation. Verification of the photon efficiencies and energy scale in data with respect to the 
simulation are discussed in Section [6] 


6 Systematic uncertainties 

Systematic uncertainties are assigned to each step of the analysis procedure using methods 
detailed in this section. Tables [3] and [4] present a summary of these uncertainties, which are 
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similar in magnitude to, or somewhat larger than the statistical uncertainties, depending on 
the photon Ej. 

The muon efficiency, taken from simulation, is corrected as a function of pj and rj using a 
method derived from data and described in Ref. m. The statistical uncertainties for these 
corrections constitute a systematic uncertainty, which we also take from Ref. Ifl7l . In addition, 
we assign a 0.5% uncertainty to account for the modifications of the standard isolation variable. 
We propagate the uncertainty in the muon efficiency by shifting the per-event weights up and 
down by one unit of systematic uncertainty. 

The photon Ej scale is potentially an important source of systematic uncertainty although 
simulations indicate that the bias in PF photon energy is negligibly small. We verify the fi¬ 
delity of the simulation by introducing an extra requirement, 0.05 < AR /<7 < 0.9, which 
gives us a high-purity subset of signal events in which the energy of the photon can be es¬ 
timated from just the muon kinematics. We refer to this estimate as E k ' n . The quantity s = 
1 — (M^ <7 — M^J/(M| — « 1 — E 7 F /E 7 m is distributed as a skewed Gaussian with a 

mean close to zero. We conducted detailed quantitative studies of the s distribution in bins of 
E FF , separately in the barrel and endcaps. We fit the distributions to a Gaussian-like function 
in which the width parameter is itself a function of s, namely, <r(s) = c(l + e bs ), with b and 
c as free parameters. Examples are given in Fig. [2] Overall, the simulation reproduces the s 
distributions in data very well. We derive some small corrections from the differences in data 
and simulation as a function of E FF and construct an alternate response matrix. The unfolded 
spectrum we obtained with this alternate response matrix differs from the original by less than 
0.2% for Ej < 40GeV, by less than 1% for E T < 75GeV, and by less than 3% in the highest 
Ej bin. We assign respective systematic uncertainties of 0.5%, 1%, and 3% for these three Ej 
ranges to account for the photon energy scale uncertainty. 


4.7 fb 1 (7 TeV) 


4.7 fb 1 (7 TeV) 



s 


s 


Figure 2: Two examples of an s distribution s = 1 — (M ^, 7 — Mj lfl )/ (M| — fit with 

a skewed Gaussian as described in the text. The left and right plots pertain to photons in 
the ECAF barrel with 5 < Ej < 10 GeV and in the ECAF endcaps with 20 < Ej < 40 GeV, 
respectively. The circle points and solid curve represent the data and the triangle points and 
dotted curve represent the simulation. 


The photon energy resolution uncertainty is well constrained by studies with electrons and 
FSR photons Q. To assess the impact of the uncertainty in the resolution, we degrade the pho¬ 
ton energy resolution in simulated events by adding in quadrature a 1 % term to the nominal 
resolution and construct a new response matrix. The differences in the unfolded spectrum rel¬ 
ative to the default response matrix are small, and we take these differences as the systematic 
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6 Systematic uncertainties 


Table 3: Relative systematic uncertainties for da /dEj (in percent). 


Kinematic 

requirement 

[GeV] 

Background 

estimation 

Muon 

efficiency 

Photon 

Ex-scale 

Photon Ex 
resolution 

Photon 

efficiency 

Pileup 

photons 

Unfolding 

Total 

0.05 < A R in < 3 

5 < E x < 10 

2.7 

3.0 

0.5 

1.0 

2.0 

1.5 

1.4 

5.1 

10 < E x < 15 

1.3 

2.5 

0.5 

0.5 

1.0 

0.4 

1.4 

3.4 

15 < E x < 20 

0.9 

2.5 

0.5 

0.5 

1.3 

0.1 

1.4 

3.3 

20 < E x < 25 

0.8 

2.7 

0.5 

0.5 

1.4 

<0.1 

1.4 

3.5 

25 < E x < 30 

0.7 

3.3 

0.5 

0.5 

1.5 

<0.1 

1.4 

4.0 

30 < Ej < 40 

1.0 

4.3 

0.5 

0.5 

1.1 

0.1 

1.4 

4.8 

40 < E x < 50 

2.9 

4.4 

1.0 

0.5 

2.8 

0.5 

1.4 

6.3 

50 < E x < 75 

7.2 

4.5 

1.0 

0.5 

2.0 

0.6 

1.4 

8.9 

75 < E x < 100 

15.3 

4.5 

3.0 

1.0 

6.9 

1.1 

1.4 

17.8 

0.05 < AR/iy < 0.5 

5 < E x < 10 

0.8 

2.1 

0.5 

1.0 

2.0 

0.1 

1.4 

3.5 

10 < E x < 15 

0.4 

2.0 

0.5 

0.5 

1.0 

<0.1 

1.4 

2.8 

15 < E x < 20 

0.3 

2.2 

0.5 

0.5 

1.3 

<0.1 

1.4 

3.1 

20 < E x < 25 

0.3 

2.5 

0.5 

0.5 

1.4 

<0.1 

1.4 

3.3 

25 < E x < 30 

0.2 

3.2 

0.5 

0.5 

1.5 

<0.1 

1.4 

3.9 

30 < E x < 40 

0.3 

4.3 

0.5 

0.5 

1.1 

<0.1 

1.4 

4.7 

40 < Ej < 50 

0.9 

3.9 

1.0 

0.5 

2.8 

<0.1 

1.4 

5.2 

50 < E x < 75 

2.3 

3.0 

1.0 

0.5 

2.0 

<0.1 

1.4 

4.6 

75 < E x < 100 

4.9 

3.1 

3.0 

1.0 

6.9 

0.8 

1.4 

9.7 


0.5 < A R n < 3 


5 < E x < 10 

6.4 

4.7 

0.5 

1.0 

2.0 

3.8 

1.4 

9.2 

10 < E T < 15 

2.8 

3.2 

0.5 

0.5 

1.0 

0.8 

1.4 

4.7 

15 < E t < 20 

1.9 

2.8 

0.5 

0.5 

1.3 

0.3 

1.4 

4.0 

20 < E x < 25 

1.7 

3.0 

0.5 

0.5 

1.4 

<0.1 

1.4 

4.0 

25 < E x < 30 

1.6 

3.4 

0.5 

0.5 

1.5 

0.1 

1.4 

4.3 

30 < E x < 40 

2.3 

4.4 

0.5 

0.5 

1.1 

0.2 

1.4 

5.3 

40 < Ex < 50 

6.5 

5.1 

1.0 

0.5 

2.8 

1.3 

1.4 

9.0 

50 < E x < 75 

16.1 

8.1 

1.0 

0.5 

2.0 

2.0 

1.4 

18.4 

75 < E x < 100 

34.5 

6.2 

3.0 

1.0 

6.9 

3.5 

1.4 

36.0 



0.05 < A R fiy 

< 3 and qj < 10 GeV 




5 < E x < 10 

1.4 

2.2 

0.5 

1.0 

2.0 

1.0 

1.4 

3.9 

10 < E x < 15 

0.6 

1.9 

0.5 

0.5 

1.0 

0.1 

1.4 

2.8 

15 < E T < 20 

0.4 

2.1 

0.5 

0.5 

1.3 

<0.1 

1.4 

3.0 

20 < E x < 25 

0.4 

2.4 

0.5 

0.5 

1.4 

<0.1 

1.4 

3.3 

25 < E x < 30 

0.5 

3.5 

0.5 

0.5 

1.5 

<0.1 

1.4 

4.1 

30 < Ex < 40 

0.6 

5.1 

0.5 

0.5 

1.1 

<0.1 

1.4 

5.5 

40 < E x < 50 

7.3 

4.7 

1.0 

0.5 

2.8 

1.0 

1.4 

9.4 

50 < E t < 75 

18.2 

8.5 

1.0 

0.5 

2.0 

4.4 

1.4 

20.8 

75 < E t < 100 

38.9 

6.4 

3.0 

1.0 

6.9 

<0.1 

1.4 

40.2 



0.05 < A Rf, y 

< 3 and qj > 50 GeV 




5 < E T < 10 

5.7 

4.2 

0.5 

1.0 

2.0 

1.8 

1.4 

7.8 

10 < E t < 15 

3.0 

3.3 

0.5 

0.5 

1.0 

0.4 

1.4 

4.9 

15 < E x < 20 

3.0 

2.8 

0.5 

0.5 

1.3 

0.3 

1.4 

4.6 

20 < E x < 25 

2.3 

2.7 

0.5 

0.5 

1.4 

0.2 

1.4 

4.2 

25 < E t < 30 

1.9 

2.6 

0.5 

0.5 

1.5 

0.1 

1.4 

3.9 

30 < Ex < 40 

2.9 

2.9 

0.5 

0.5 

1.1 

0.3 

1.4 

4.6 

40 < E t < 50 

1.5 

2.8 

1.0 

0.5 

2.8 

0.2 

1.4 

4.6 

50 < E x < 75 

3.9 

2.8 

1.0 

0.5 

2.0 

0.3 

1.4 

5.5 

75 < E x < 100 

8.2 

3.5 

3.0 

1.0 

6.9 

0.2 

1.4 

11.8 
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Table 4: Relative systematic uncertainties for da /d AK )r) . (in percent). 


Kinematic 

requirement 

Background 

estimation 

Muon 

efficiency 

Photon 

Er-scale 

Photon Et 
resolution 

Photon 

efficiency 

Pileup 

photons 

Unfolding 

Total 

E t > 5.0 GeV 

0.15 < A R n < 0.1 

0.7 

2.4 

<0.1 

<0.1 

1.0 

<0.1 

1.4 

3.0 

0.1 < A R,, 7 < 0.15 

0.6 

2.3 

<0.1 

<0.1 

1.1 

<0.1 

1.4 

3.0 

0.15 < AR,, 7 < 0.3 

0.4 

2.3 

<0.1 

<0.1 

1.0 

<0.1 

1.4 

2.9 

0.3 < A R,, 7 < 0.5 

0.5 

2.3 

<0.1 

<0.1 

1.0 

0.1 

1.4 

3.0 

0.5 < A R )l7 < 0.8 

1.1 

2.6 

<0.1 

<0.1 

1.0 

0.6 

1.4 

3.4 

0.8 < A R, l7 < 1.2 

2.2 

3.2 

<0.1 

<0.1 

1.1 

1.1 

1.4 

4.4 

1.2 < A R /l7 < 1.6 

4.1 

3.7 

<0.1 

<0.1 

1.1 

1.7 

1.4 

6.1 

1.6 < A R,, 7 < 2.0 

6.6 

4.9 

<0.1 

<0.1 

1.1 

2.8 

1.4 

8.8 

2.0 < AR,, 7 < 3.0 

18.3 

9.9 

<0.1 

<0.1 

1.3 

7.9 

1.4 

22.3 

£ t > 5.0 GeV and q T < 10 GeV 

0.15 < A R n < 0.1 

0.2 

2.1 

<0.1 

<0.1 

1.0 

<0.1 

1.4 

2.8 

0.1 < AR ll7 < 0.15 

0.2 

2.2 

<0.1 

<0.1 

1.1 

<0.1 

1.4 

2.8 

0.15 < AR,, 7 < 0.3 

0.1 

2.1 

<0.1 

<0.1 

1.0 

<0.1 

1.4 

2.7 

0.3 < A R, l7 < 0.5 

0.3 

2.2 

<0.1 

<0.1 

1.0 

0.1 

1.4 

2.8 

0.5 < A R t , 7 < 0.8 

0.7 

2.4 

<0.1 

<0.1 

1.0 

0.3 

1.4 

3.0 

0.8 < A R, l7 < 1.2 

1.3 

2.5 

<0.1 

<0.1 

1.1 

0.6 

1.4 

3.4 

1.2 < AR fl7 < 1.6 

2.2 

2.7 

<0.1 

<0.1 

1.1 

1.0 

1.4 

4.1 

1.6 < A R n < 2.0 

3.8 

3.1 

<0.1 

<0.1 

1.1 

2.1 

1.4 

5.6 

2.0 < A R,, 7 < 3.0 

15.9 

7.4 

<0.1 

<0.1 

1.3 

9.0 

1.4 

19.8 

£ t > 5.0 GeV and q T > 50 GeV 

0.15 < ARj, 7 < 0.1 

1.8 

2.5 

<0.1 

<0.1 

1.0 

<0.1 

1.4 

3.6 

0.1 < A R,, 7 < 0.15 

1.1 

2.3 

<0.1 

<0.1 

1.1 

<0.1 

1.4 

3.1 

0.15 < A R,, 7 < 0.3 

1.5 

2.4 

<0.1 

<0.1 

1.0 

<0.1 

1.4 

3.4 

0.3 < A R, l7 < 0.5 

1.7 

2.4 

<0.1 

<0.1 

1.0 

0.1 

1.4 

3.4 

0.5 < A R ll7 < 0.8 

2.6 

2.9 

<0.1 

<0.1 

1.0 

0.7 

1.4 

4.4 

0.8 < AR ll7 < 1.2 

4.2 

3.8 

<0.1 

<0.1 

1.1 

1.4 

1.4 

6.1 

1.2 < AR,, 7 < 1.6 

9.1 

5.2 

<0.1 

<0.1 

1.1 

1.9 

1.4 

10.8 

1.6 < A R n < 2.0 

14.9 

7.4 

<0.1 

<0.1 

1.1 

3.4 

1.4 

17.1 

2.0 < A R,, 7 < 3.0 

22.3 

10.3 

<0.1 

<0.1 

1.3 

5.1 

1.4 

25.1 
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7 Results 


uncertainty due to photon energy resolution. 

Efficiency corrections for photons are applied as part of the unfolding procedure described in 
Section [5] and are derived from the simulation. We verify these corrections using the data in 
the following way. An isolated FSR photon with Ej > 5 GeV nearly always produces a cluster 
in the ECAL. We define an efficiency to reconstruct and select PF photons given such isolated 
clusters. This efficiency rises from 60% for Ej between 5-10 GeV to approximately 90% for 
Ej > 50 GeV and is nearly the same in data and simulation. We take the difference, added in 
quadrature to the statistical uncertainties of the efficiencies, as the systematic uncertainty. 

As described briefly in Section [5J the unfolding procedure has been cross-checked in several 
ways. To assess a systematic uncertainty due to unfolding, we use the small discrepancies 
observed in the closure test. 

The uncertainty in the background estimate is dominated by the uncertainties associated with 
the corrections that we obtained from the control region (Section |4jl. The statistical uncertainty 
in the weights for jet multiplicity has a negligible impact, as does the correction for charged 
hadrons in the photon isolation cone. The parameterized functions to correct the photon dis¬ 
tributions in Ej and // carry statistical uncertainties that we propagate to the measured cross 
sections through simplified MC models. Since the nonprompt photon Ej, //, and A/% 7 distribu¬ 
tions in the control and signal regions are indistinguishable, we do not assess any uncertainty 
in the modeling. 

The uncertainties in the non-DY backgrounds (tt and diboson production) are obtained from 
the uncertainties in the theoretical cross sections, the luminosity, and the statistical uncertainty 
in the simulated event samples. We assign 50% uncertainty to the W+jets and multijet back¬ 
ground estimates, which are quite small. 

The systematic uncertainty from the simulation of pileup depends primarily on the assumed 
cross section for additional pp collisions (roughly the same as the minimum-bias cross sec¬ 
tion) 1251 . We vary the value of this cross section by 5% and evaluate the impact on the unfolded 
spectra. 

The uncertainty in the integrated luminosity is 2.2% li26l . 

Theoretical uncertainties have been calculated and pertain to the reported theoretical predic¬ 
tion only. We propagated the uncertainty due to parton distribution functions (PDFs) using 
the prescription of Ref. Il27l . We vary the factorization/renormalization scale parameters by a 
factor of 2 to estimate associated scale uncertainties introduced due to neglected higher-order 
quantum corrections. Finally, we include the MC statistical uncertainty. 


7 Results 

The differential cross sections are obtained by subtracting the estimated backgrounds from the 
observed distributions, unfolding the result, and dividing by the bin width and the integrated 
luminosity, C = 4.7 fb 1 . No acceptance correction is applied, so these cross sections are de¬ 
fined relative to the kinematic and fiducial requirements listed in Table [T] 

The measured differential cross sections da/d Ej and da / dAR, n are displayed in Fig. |3] and 
listed in Tables [5] and [6] A bin-centering correction is applied following the method of Reh H28l ; 
the abscissa of each point is based on the integral of the simulation across the bin and on the bin 
width. The shaded region represents the prediction and uncertainty from POWHEG+PYTHIA, 
obtained at the parton level: only the requirements in Table[l]have been applied to the generator- 
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level muons and photons. The agreement with data is good. 

Energy spectra for photons closer to (0.05 < A< 0.5) and farther from the muon (0.5 < 
A Rp l7 < 3) are shown in Fig. [4] The rates for photons with large Aand Ej are also well 
reproduced. The number of events with 30 < < 87 GeV is about 18% of the number with 

60 < M w< < 120 GeV. Of the events with 30 < M ?/?( < 87 GeV, the fraction of events with at 
least one photon with Ej > 5 GeV and 0.05 < AR^ 7 < 0.5 is 8.7 ± 0.1 (stat) ± 0.2 (syst)%, and 
with 0.5 < A R tl y < 3 is 5.6 ± 0.1 (stat) ± 0.2 (syst)%. Photons with AR^ > 1.2 and Ej > 40 GeV 
constitute a small fraction (1.3 ± 0.5 (stat) ± 0.6 (syst)) x 10~ 4 . 


4.7 fb' 1 (7 TeV) 



E t [GeV] 



4 

| 2 
Q 0 
2 -2 
m -A 


• 

•• • • * 

: 

0 0.5 

1.5 2 2.5 3 


AR^ 


Figure 3: Measured differential cross sections dcr/dEj (left) and dcr/dAR^ (right). In the 
upper panels, the dots with error bars represent the data, and the shaded bands represent the 
POWHEG+PYTHIA calculation including theoretical uncertainties. The central panels display 
the ratio of data to the MC expectation. The lower panels show the standard deviations of the 
measurements with respect to the calculation. A bin-centering procedure has been applied. 

We define two subsamples of signal events, one with the Z boson transverse momentum qj < 
10 GeV, and the other with qj > 50 GeV. The measured cross sections shown in Fig. [^demon¬ 
strate rather different energy spectra for these two cases, though dcr /dAis basically the 
same. 

As a final illustration of the nature of this event sample, we present distributions of dimuon 
mass (Mmi) and the three-body mass (M mn ) in Fig. [ 6 j The small increase in the ratio of data to 
theory for < 40 GeV reflects the insufficient next-to-leading-order accuracy of the simula¬ 
tion; the kinematic requirements on the muons induce a loss of acceptance that require higher- 
order QCD corrections, as discussed in Ref. Ifl7| . Although the masses of the dimuon pairs 
populate the tail of the Z resonance (in fact they were selected this way), the three-body mass 
distribution displays a nearly-symmetric resonance peak at the mass of the Z boson, thereby 
confirming the identity of these events as radiative decays Z —» u 1 fi ~ 7 . 
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7 Results 


Table 5: Measured differential cross section da/ d£ T in pb/GeV. For the data values, the first 
uncertainty is statistical and the second is systematic. For the theory values, the uncertainty 
combines statistical, PDF, and renormalization/factorization scale components. 

Kinematic requirement [GeV] 

Data 

POWHEG+PYTHIA 

0.05 < A R, n < 3 

5 < E T < 10 


1.260 ±0.015 ±0.070 

1.270 ±0.075 

10 < E T < 15 


0.685 ± 0.009 ± 0.028 

0.694 ± 0.040 

15 < E T < 20 


0.411 ± 0.006 ± 0.016 

0.433 ± 0.025 

20 < E T < 25 


0.267 ±0.005 ±0.011 

0.280 ±0.017 

25 < E t < 30 


0.170 ± 0.004 ± 0.008 

0.177 ± 0.011 

30 < Ex < 40 


(7.26 ± 0.19 ± 0.39) x 10^ 2 

(7.20 ± 0.42) x 10~ 2 

40 < Ej < 50 


(1.49 ± 0.09 ± 0.10) x 10^ 2 

(1.34 ± 0.08) xl0~ 2 

50 < E t < 75 


(2.68 ± 0.26 ± 0.25) x lO^ 3 

(2.27 ± 0.14) x 10“ 3 

75 < E T < 100 


(5.81 ± 1.16 ± 1.00) x 10^ 4 

(3.47 ± 0.32) x 10~ 4 

0.05 < A R flJ < 0.5 

5 < E t < 10 


0.749 ± 0.009 ± 0.031 

0.779 ± 0.045 

10 < E t < 15 


0.417 ±0.006 ±0.015 

0.433 ± 0.025 

15 < E t < 20 


0.256 ± 0.005 ± 0.010 

0.272 ±0.016 

20 < E t < 25 


0.168 ± 0.004 ± 0.007 

0.177 ± 0.011 

25 < E t < 30 


0.105 ± 0.003 ± 0.005 

0.112 ±0.007 

30 < Ex < 40 


(4.51 ± 0.14 ± 0.23) x 10~ 2 

(4.44 ± 0.26) x 10~ 2 

40 < Ex < 50 


(8.93 ± 0.65 ± 0.51) x lO^ 3 

(8.53 ± 0.50) xl0~ 3 

50 < E t < 75 


(1.80 ± 0.18 ± 0.09) x 10- 3 

(1.63 ± 0.10) xl0~ 3 

75 < E t < 100 


(3.58 ± 0.98 ± 0.36) x 10^ 4 

(2.42 ± 0.37) x 10~ 4 

0.5 < A R in < 3 

5 < E t < 10 


0.513 ± 0.012 ± 0.049 

0.489 ± 0.028 

10 < E t < 15 


0.268 ± 0.006 ± 0.014 

0.260 ±0.015 

15 < E t < 20 


0.155 ± 0.004 ± 0.007 

0.161 ±0.010 

20 < E t < 25 


(9.94 ± 0.33 ± 0.45) x 10~ 2 

(1.03 ± 0.06) xlO- 1 

25 < E t < 30 


(6.52 ± 0.26 ± 0.32) x lO^ 2 

(6.55 ± 0.39) x 10 -2 

30 < Ex 5 40 


(2.76 ± 0.12 ± 0.16) x 10 -2 

(2.76 ± 0.17) x 10 -2 

40 < Ex < 50 


(6.01 ± 0.67 ± 0.56) x 10 -3 

(4.85 ± 0.33) x 10~ 3 

50 < E t < 75 


(8.75 ± 1.86 ± 1.60) x 10 -4 

(6.38 ± 0.60) xl0~ 4 

75 < E t < 100 


(2.23 ± 0.63 ± 0.80) x 10^ 4 

(1.04 ± 0.27) x 10“ 4 

0.05 < A R fn < 3 and qj < 10 GeV 

5 < E t < 10 


0.527 ±0.009 ±0.024 

0.535 ± 0.033 

10 < E t < 15 


0.294 ± 0.005 ± 0.010 

0.296 ±0.018 

15 < E t < 20 


0.184 ± 0.004 ± 0.007 

0.191 ±0.012 

20 < E t < 25 


0.127 ±0.003 ±0.005 

0.129 ± 0.008 

25 < E t < 30 


(8.59 ± 0.28 ± 0.40) x 10^ 2 

(8.25 ± 0.54) x 10~ 2 

30 < Ex < 40 


(3.22 ± 0.12 ± 0.19) x 10 -2 

(2.89 ± 0.18) x 10~ 2 

40 < E T < 50 


(1.46 ± 0.27 ± 0.14) x 10- 3 

(1.14 ± 0.12) x 10~ 3 

50 < E t < 75 


(1.92 ± 0.67 ± 0.42) x 10- 4 

(8.44 ± 1.60) x 10~ 5 

75 < E t < 100 


(1.67 ± 2.10 ± 0.66) x 10^ 5 

(6.66 ± 5.13) x 10~ 6 

0.05 < A R tn < 3 and qj > 50 GeV 

5 < E t < 10 


0.104 ± 0.005 ± 0.008 

0.095 ± 0.005 

10 < E t < 15 


(6.26 ± 0.28 ± 0.33) x lO^ 2 

(5.72 ± 0.31) x 10~ 2 

15 < E t < 20 


(3.67 ± 0.20 ± 0.19) x 10 -2 

(3.38 ± 0.18) xl0~ 2 

20 < E t < 25 


(2.19 ± 0.15 ± 0.10) x 10^ 2 

(2.32 ± 0.13) x 10~ 2 

25 < E t < 30 


(1.94 ± 0.14 ± 0.09) x 10- 2 

(1.64 ± 0.10) x 10~ 2 

30 < Ex T 40 


(9.98 ± 0.71 ± 0.51) x 10^ 3 

(9.79 ± 0.55) xl0~ 3 

40 < Ex < 50 


(6.21 ± 0.55 ± 0.32) x 10~ 3 

(5.58 ± 0.33) x 10~ 3 

50 < E t < 75 


(1.90 ± 0.20 ± 0.11) x 10^ 3 

(1.76 ± 0.11) x 10“ 3 

75 < E t < 100 


(4.56 ± 0.95 ± 0.55) x 10~ 4 

(3.13 ± 0.30) xl0~ 4 
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Table 6 : Measured differential cross section da /d AR „ 7 in pb. For the data values, the first 
uncertainty is statistical and the second is systematic. For the theory values, the uncertainty 
combines statistical, PDF, and renormalization/factorization scale components. 


Kinematic 

requirement 




Data 

POWHEG+PYTHIA 

Ft > 5.0 GeV 

0.05 

< 

AR^ij 

< 

0.1 

53.90 

±0.76 

± 2.00 

56.60 

± 

3.26 

0.1 

< 

A 

< 

0.15 

31.90 

±0.59 

± 1.20 

33.20 

± 

1.96 

0.15 

< 

A 

< 

0.3 

18.40 

±0.25 

±0.67 

19.00 

± 

1.10 

0.3 

< 

A Rjij 

< 

0.5 

10.10 

±0.16 

±0.37 

10.50 

± 

0.59 

0.5 

< 

A Rjiy 

< 

0.8 

6.14 

± 0.11 

±0.25 

6.29 

± 

0.37 

0.8 

< 

A Rjiy 

< 

1.2 

4.22 

±0.09 

± 0.21 

4.10 

± 

0.24 

1.2 

< 

A Rjiy 

< 

1.6 

2.94 

±0.08 

±0.19 

2.91 

± 

0.17 

1.6 

< 

A Rjij 

< 

2.0 

1.76 

±0.07 

±0.16 

1.79 

± 

0.11 

2.0 

< 

ARji'j- 

< 

3.0 

0.46 

±0.04 

± 0.10 

0.33 

± 

0.02 






Ej > 

GeV and qj < 10 GeV 



0.05 

< 

ar ,, 7 

< 

0.1 

23.00 

±0.50 

±0.82 

24.40 

± 

1.53 

0.1 

< 

ar ,, 7 

< 

0.15 

13.70 

±0.39 

±0.49 

14.20 

± 

0.88 

0.15 

< 

ar ,, 7 

< 

0.3 

7.88 

±0.17 

±0.28 

8.21 

± 

0.51 

0.3 

< 

ar ,, 7 

< 

0.5 

4.38 

± 0.10 

±0.16 

4.48 

± 

0.28 

0.5 

< 

ar ,, 7 

< 

0.8 

2.65 

±0.07 

± 0.10 

2.67 

± 

0.17 

0.8 

< 

ar ,, 7 

< 

1.2 

1.75 

±0.05 

±0.07 

1.75 

± 

0.11 

1.2 

< 

ar ,, 7 

< 

1.6 

1.29 

±0.05 

±0.06 

1.25 

± 

0.08 

1.6 

< 

ar ,, 7 

< 

2.0 

0.72 

±0.04 

±0.04 

0.79 

± 

0.05 

2.0 

< 

ar,, 7 

< 

3.0 

0.10 

± 0.01 

± 0.02 

0.09 

± 

0.01 






E T > 5 

.OGeV 

and qj > 50 GeV 



0.05 

< 

A R ;< 7 


0.1 

4.94 

±0.23 

± 0.21 

5.07 

± 

0.27 

0.1 

< 

A R ;< 7 

< 

0.15 

2.97 

±0.18 

± 0.11 

3.05 

± 

0.18 

0.15 

< 

A R ;< 7 

< 

0.3 

1.71 

±0.08 

±0.07 

1.73 

± 

0.09 

0.3 

< 

ar ,, 7 

< 

0.5 

0.95 

±0.05 

±0.04 

0.98 

± 

0.06 

0.5 

< 

AR|i 7 

< 

0.8 

0.62 

±0.04 

±0.03 

0.58 

± 

0.03 

0.8 

< 

AR|i 7 

< 

1.2 

0.44 

±0.03 

±0.03 

0.37 

± 

0.02 

1.2 

< 

AR|i 7 

< 

1.6 

0.22 

±0.03 

± 0.02 

0.19 

± 

0.01 

1.6 

< 

ARji 7 

< 

2.0 

0.13 

± 0.02 

± 0.02 

0.10 

± 

0.01 

2.0 

< 

ar ,, 7 

< 

3.0 

(8.45 

±1.38 

± 2 . 10 ) x 10 

~ 2 (3.62 

± 

0.24) x 10 ^ 2 
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8 Summary 


4.7 fb' 1 (7 TeV) 


4.7 fb' 1 (7 TeV) 




4 

S 2 
Q n 

• • 

• 

• 


4 

S 2 
Q 0 

• 

• 

• 

• • • 

• • 

2 -2 
w -4 





2 -2 
w -4 






10 20 30 40 50 60 70 80 9 

0 100 

10 2 

0 30 40 50 60 70 80 90 1C 


E t [GeV] 


E t [GeV] 


Figure 4: Measured differential cross sections da/dEj for photons close to the muon (0.05 < 
AR^ 7 < 0.5, left) and far from the muon (0.5 < A < 3, right). The dots with error bars 
represent the data, and the shaded bands represent the POWHEG+PYTHIA calculation including 
theoretical uncertainties. The central panels display the ratio of data to the MC expectation. The 
lower panels show the standard deviations of the measurements with respect to the calculation. 
A bin-centering procedure has been applied. 

8 Summary 

A study of final-state radiation in Z boson decays was presented. This study serves to test the 
simulation of events where mixed QED and QCD corrections are important. The analysis was 
performed on a sample of pp collision data at \Js = 7 TeV recorded in 2011 with the CMS de¬ 
tector and corresponding to an integrated luminosity of 4.7 fb -1 . Events with two oppositely 
charged muons and an energetic, isolated photon were selected with only modest backgrounds. 
The differential cross sections d a /d Ej and da /d AE irj , were measured for photons within the 
fiducial and kinematic requirements specified in Table [lj and comparisons of da /dEj for pho¬ 
tons close to a muon and far from both muons were made. In addition, the differential cross 
sections da/dEj and da/ dAR f , 7 were compared for events with large and small transverse 
momentum of the Z boson, as computed from the two muons and the photon. Simulations 
based on POWHEG+PYTHIA reproduce the CMS data well, with discrepancies below 5% for 
5 < Ej < 50 GeV and 0.05 < AR ? , 7 < 2 as quantified in Tables [5] and [6j 
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Figure 5: Measured differential cross sections da/dEj and da /d AK ;rj , for qj < 10 GeV (top 
row) and qj > 50 GeV (bottom row). The dots with error bars represent the data, and the 
shaded bands represent the POWHEG+PYTHIA calculation including theoretical uncertainties. 
The central panels display the ratio of data to the MC expectation. The lower panels show 
the standard deviations of the measurements with respect to the calculation. A bin-centering 
procedure has been applied. 
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tation. The lower panels show the standard deviations of the measurements with respect to the 
calculation. A bin-centering procedure has been applied. 
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